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INTRODUCTION 

Systems which are overdamped due to the frictional drag of solvent 
have the motion of their constituents described by the diffusion equation 
as an approximation to the Fokker-Planck equation (!)• Because of the 
frequent occurrence of this condition, diffusional models are commonly 
used to describe biological, physical and chemical processes including 
heat conduction, mass transport, some aspects of chemical reactions, 
and so on. In order to apply the general diffusion equation to a particular 
physical situation or model, the space in which diffusion occurs must be 
designated and the behavior of appropriate combinations of the particle 
density and its derivatives specified at the boundaries. However, in the 
physical situation in which there is partial or complete absorption 
(capture or reaction depending on the particular problem being modeled) 
at the boundary, it has not been completely clear how to interpret the 
parameters determining the boundary condition. In Ref. 2 in the context 
of the diffusion-collision model of protein folding, two possible 
parameter choices were discussed in detail and a third possibility, 
the one to be outlined below in the present paper, was mentioned without 
proof. Since the third possibility leads to the rates of diffusion- 
controlled processes having a solvent viscosity-independent part as 
well as a viscosity-dependent part, it has important consequences 
for intramolecular and intermolecular movements in biological systems. 

In addition, recent numerical studies (3) and approximate analytical 
calculations (4) on one-dimensional diffusional systems have lent 
credence to the third possibility. Therefore, the purpose of the 
present paper is to present a derivation of the "radiation” or partially 
absorbing boundary condition (the perfect absorber being a special 
case) for a diffusion-controlled process in which not every collision 

W85>.- 
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between a Brownian particle and the absorbing boundary results in 
removal of the particle from the diffusion space. The details of 
the derivation are follwed by a discussion of some experiments Involving 
Intramolecular polypeptide movements in which the solvent viscosity 
was varied, and the implications for diffusion in membranes in which 
localization due to trapping takes place (5,6). 


THE DERIVATION 

To be specific, one-dimensional diffusion will be discussed below. 

The designation of the space determines whether the system is infinite 
or finite, and in the latter case one must specify the length of 
the space in order to calculate quantities of physical interest in 
absolute terms. 

Consider finite space diffusion first (infinite space diffusion 
will be considered later). Boundary conditions must be specified 
at two boundaries, x = a and x = b with b > a. Let the particle density 
be p(x,t) at position x and time t. The diffusion equation is 

H ^ a < X < b (1) 

with D the diffusion coefficient (assumed constant). The boundary 
conditions at a and b, which along with the initial value p(x,o) are 
needed to specify p(x,t), are to be linear combinations of p and its 
first derivatives Sp/3x and 3p/3t. In the absence of particle emission 
at the boundaries, 3p/3t will not contribute to the boundary conditions. 

A 

Since the particle flux is ^(x,t) = -D3p/3x x, which is the net number 
of particles per second crossing a unit area perpendicular to the x-axis. 
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a completely reflecting boundary at which no absorption takes place is 
uniquely specified by the condition ^ = 0. Let this boundary condition 
prevail at x » b in the one-dimensional finite space example. Thus, 

3p/3x = 0 at X » b. 

At the boundary x = a, let the conditions be such that partial 
absorption takes place. That is, some fraction of the Brownian particles 
that reach the boundary at x = a are not reflected and no longer diffuse. 
This general case encompasses the limits of perfect absorption and 
perfect reflection as well as all intermediate situations. For partial 
absorption, the reflected flux at x » a will be smaller than the Incident 
flux by the number of particles per second being absorbed, that is 
3^ ■ y3_, 0 ^ Y £ 1, where 3_ is the magnitude of the flux toward 
the boundary at x » a and 3^ is the magnitude of the reflected flux 
at X = a (that is, the portion of the net flux at x = a going in the 
positive X direction). In the spirit of kinetic theory, one may write 
3^(a,t) = <v>^p^(a,t) where p^(a,t) are the reflected (+) and incident 
(-) particle densities and <v>^ are "average" values for the reflected 
and incident particle velocities at x = a derived from the mean value 
theorem of the calculus. If the absorption mechanism at the boundary 
is Independent of velocity, then <v>^ = <v> = <v>. This will be 

assumed as the simplest possibility. In many situations the boundary 
will be due to a potential energy barrier whose height sets a minimum 
value for the kinetic energy of the incident particles in order to 
be absorbed at x = a. Nevertheless, if the approach to a Maxwell- 
Boltzmann velocity distribution is fast on the diffusion time scale. 
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then the incident and reflected particle velocity distributions will be 
virtually identical. Then, since 

3(a,t) = 3_|_(a,t) - 3_(a,t) = (y - l)<v>p_(a,t) (2) 

and 


p(a,t) « p^(a,t) + p_(a,t) = (1 + y)P_(a,t) 

one obtains the boundary condition 
3(a,t) 
or 


~a V y) 


3p 

3x- 


(1 - y) <v> 
(1 + y) D 


(3) 


(4) 


(5) 


Note that with this boundary condition, perfect absorption y = 0 does not 

imply p| “0. However, for perfect reflection y = 1, the condition 
a 

3p/3x| = 0 is found* The consequence of Eq* (5) for several biophysical 

a 

situations have been briefly mentioned by the present author previously 
(2,6,7) without discussion# Below some of the implications and relevant 
experiments are analyzed# 


THE IMPLICATIONS 


Eq# (5) has important implications about the rates of some diffusional 
processes# The rates of absorption in the model outlined above may be 
calculated using an extension of the idea of first passage times (8), as 
shown by a number of authors (9,10,11). This is carried out below for 
the diffusional model considered above by writing the fraction of particles 
absorbed at x = a in terms of the mean absorption time (Inverse of the 
absorption rate. 
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Let the fraction of absorbed particles be n (t), so that 

* 

D^(t) - 1 - J dx p(x,t) (6) 

where p Is normalized Initially to one particle (thus, p may be thought 
of as the relative position probability density). Solution of the 
diffusion equation with the given boundary conditions and an Initial 
distribution In which the particles start at x “ b Indicates that the 
Infinite series of decreasing exponential terms which represents the 
time dependence of n (t) may be well-approximated by a single exponential. 
That Is, 



Y ■ 0 perfect absorber situation does not (using Eq. (5)) lead to 
p| ■ 0 as the boundary condition. It does, however, lead to the boundary 
condition p(x,t) *0atx*a-A»A>0 where A is a length parameter 
(referred to as the Milne extrapolation length) whose value is 
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A . D 

^ (1-Y) <v> 


(9) 


Eq. (9J was found by computing the rate of particle injection at x * b 
needed to have the above diffusional model at steady-state with one particle 
in the region (a,b) in order to compare with the results of Refs. 3 and A. 
This comparison also allows one to more precisely identify <v> to be 
(a/2)^ where o = k„T/m with k„ the Boltzmann constant, T the absolute 

"-is JJ 

temperature and m the Brownian particle mass. Note that a is the mean 
square velocity in one-dimension for a Maxwell-Boltzmann distribution. 
Another interpretation is <v> = where is the most probable 

speed. 


SOME RELEVANT EXPERIMENTS 

The "radiation” boundary condition derived above leads to the possibility 
that the rates of diffusion driven processes may not be directly pro- 
portional to D and, in fact, may be effectively independent of D for 
small enough y* this regard, some experiments have been carried out 

to study the intramolecular motion of polypeptide chain segments in 
solvents of different viscosities (since, according to the Einstein 
relation D = k^T/f with f the friction coefficient, varying the viscosity 

D 

ought to vary D in an inverse manner). Haas et al. (12) investigated the 
kinetics of the fluorescence decay of the energy donor in a homologous 
series of oligopeptides, each containing at its ends a donor and an 
acceptor of electronic excitation decay energy in solvent mixtures of 
different viscosities. With an assumed theoretical analysis, diffusion 
coefficients were derived which increased systematically upon decreasing 
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the solvent viscosity. The values obtained for the diffusion coefficients 
were about an order-of-magnitude smaller than the values expected for the 
diffusion coefficients of the free chromophores in solvents of comparable 
viscosity, and appear to have a solvent viscosity independent part in 
the derived values of D Although the results obtained for D in this 
experiment are probably model dependent, there is a clear dependence 
of diffusion coefficient on solvent viscosity, as well as the suggestion 
of a viscosity independent contribution to D In the analysis of this 
experiment the boundary condition p (a,t) = 0 was used rather than 
Eq. (5) with y “ 0* 

Tsong and Baldwin (13), on the other hand, in their study of the 
kinetics of folding of the two forms of unfolded ribonuclease A (with all 
disulfide bonds intact) as a function of solvent viscosity, by adding 
either sucrose or glycerol, found no dependence on solvent viscosity, 
the rates of the two observed reactions being either unchanged or slightly 
faster in the presence of sucrose or glycerol. In the same system, 

Tsong (14) has recently found a reaction which is strongly dependent 
on solvent viscosity and somewaht faster than the reactions observed 
in Ref. 13. 

The precise interpretation of these experimental results is not 
completely clear, at present, due to the complexity of the systems 
involved. However, it appears that diffusion mediated reactions are 
playing a signfiicant role, and that there are both solvent viscosity 
dependent and solvent viscosity Independent contributions to the kinetics, 
perhaps of the kind suggested by Eq. (8) and the subsequent discussion. 

If one assumes that considerations similar to those leading to 
Eq. (5) apply also to the two-dimensional surface of membranes, then 
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localization in such systems due to diffusion-driven trapping (5,6) 
will also show the mixed dependence on D of Eq, (8), although the detailed 
form will be different (see Ref. 6). In that case, D would not be 
simply directly proportional to the localization rate. This problem 
requires further study. Some approximate results have been obtained 
by Harris (15) for the related problem of steady absorption of Brownian 
particles by an infinite right circular cylinder, corresponding to 
absorption by a circular sink in two dimensions. A modified boundary 
condition was again found by introducing an effective radius of absorption 
corresponding to the Milne extrapolation length in one-dimension. 

The extension of the above ideas to infinite systems leads to 
further interesting results concerning approach to equilibrium effects. 
Consider, as a definite example, a spherically symmetrical infinite 
system of molecules of initial uniform concentration surrounding 
a target of radius a centered at r = o. The molecules diffuse in 
the infinite space, and their concentration as a function of position 
and time, p(r,t) satisfies the diffusion equation, Eq. (1), within 
the limits a < r < Corresponding to Eq. (5), the boundary condition 
appropriate to the infinite space is 


aL ^ 
3r 


■a 

where aL will be the three-dimensional analog of (1 + Y)D/((1 - Y)<v>) 
in Eq. (5). As r -► », p(r,t) p^ for the other boundary condition. 

The above diffusion problem may be solved, in principle, using 
the Laplace transform method to obtain Aua^D |^, the time dependent 
rate of association. To proceed further, suppose that the initial 


( 10 ) 
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concentration of targets Is equal to also. Then, the theory of 
von Smoluchowskl (16) for colloidal particles may be applied here 
to yield n (t) the probability that a target molecule pair Is formed 
before time t. The result Is 


n^(t) 


(t/I(l -I- DTD 
1 + (t/[(l+L)T]) 


where T s l/(47raDp^), the "unit" 


of time appropriate to this physical 


( 11 ) 


situation. Again, there appears a mixed D dependence in the overall 
rate of association. Harris (17) has recently studied numerically 
the related problem of steady-state absorption by a sphere. 
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